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In the presence of the Josephson vortex lattice in layered superconductors, a 
small c-axis magnetic field penetrates in the form of vortex chains. In gen- 
eral, the structure of a single chain is determined by the ratio of the London 
[X] and Josephson [Xj] lengths, a = X/Xj. The chain is composed of tilted 
vortices at large a's ( tilted chain ) and at small a 's it consists of a cross- 
ing array of Josephson vortices and pancake-vortex stacks (crossing chain). 
We study chain structures at intermediate a 's and found two types of phase 
transitions. For a < 0.6 the ground state is given by the crossing chain in a 
wide range of pancake separations a > [2 — 3]Aj. However, due to attractive 
coupling between deformed pancake stacks, the equilibrium separation can 
not exceed some maximum value depending on the in-plane field and a. The 
first phase transition takes place with decreasing pancake-stack separation a 
at a = [1 — 2]Xj, and rather wide range of the ratio a, 0.4 < a ^ 0.65. With 
decreasing a, the crossing chain goes through intermediate strongly- deformed 
configurations and smoothly transforms into a tilted chain via a second-order 
phase transition. Another phase transition occurs at very small densities of 
pancake vortices, a ~ [20 — 30] A j, and only when a exceeds a certain crit- 
ical value ~ 0.5. In this case a small c-axis field penetrates in the form of 
kinks. However, at very small concentration of kinks, the kinked chains are 
replaced with strongly deformed crossing chains via a first-order phase tran- 
sition. This transition is accompanied by a very large jump in the pancake 
density. 

PACS numbers: 74.25. Qt, 74.25. Op, 74.20. De 



1. INTRODUCTION 



Layered superconductors have an amazingly rich phase diagram in tilted 
magnetic field. Point (or "pancake") vortices generated by the c-axis com- 
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Fig. 1. Crossing (left) and tilted (right) vortex chains. 



ponent of the magnetic field ^ can form a vast variety of lattice structm'es. 
Possible structures include the kinked lattice, l -^ l ^ l ^ l ^ l, tilted vortex chains 
coexisting lattices with different orientatioiP and crossing lattices composed 
of sublattices of Josephson vortices (JVs) and pancake-vortex stack^^. The 
main source of such richness is the existence of two very different kinds of 
interactions between the pancake vortices in different layers: magnetic and 
Josephson interactions. The key parameter, which determines the relative 
strength of these two interactions and plays a major role in selecting the 
lattice structures, is the ratio of the two fundamental lengths, the in-plane 
London penetration depth, A = Xab, and Josephson length Aj = 7s, with 7 
being the anisotropy parameter and s being the interlayer spacing, a = A/Aj. 
One can distinguish two limiting cases which we refer to as "extremely 
anisotropic" case, a < 0.4, and "moderately anisotropic" case a > 0.7. 
Among known atomically layered superconductors, only Bi2Sr2CaCu20i. 
(BSCCO) and related compounds may belong to the "extremely anisotropic" 
family. Even in this compound the parameter a is not smaller than ~ 0.25 
and increases with temperature so that BSCCO typically becomes "moder- 
ately anisotropic" in the vicinity of the transition temperature. 

At very small c-axis fields (up to 1-2 G) the pancake stacks in layered 
superconductors within a wide range of anisotropics are arranged in chains. 
An isolated chain is a two-dimensional array of pancake vortices oriented 
perpendicular to the layers, see Fig. ^ At somewhat higher fields the chains 
are surrounded by the regions of regular vortex lattic^. The internal struc- 
ture of an isolated chain depends on the ratio a and it is relatively simple 
in two limiting cases. At large a the chain is composed of tilted pancake 
stacks (tilted chain, right picture in Fig. ^ and at small a it consists of 
crossing arrays of JVs and pancake stacks (crossing chains, left picture in 
Fig. ^ . In this paper we address the highly nontrivial problem of how one 
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structure transforms into another in the region of intermediate a. We analyt- 
icaUy and numerically computed ground state configurations in the isolated 
vortex chain and found two types of phase transitions. The first phase tran- 
sition typically takes place for the intermediate separations between pancake 
stacks a, a = [1 — 2]Aj, and rather wide range of the ratio a, 0.4 < a < 0.65. 
The ground state is given by the crossing chain in a wide range of pancake 
separations a. Due to attraction between deformed pancake stacks the 
equilibrium separation can not exceed some maximum value, which depends 
on the in-plane field and a and it is typically of the order of several Aj. With 
decrease of the pancake separation a, the crossing chain becomes strongly 
deformed and smoothly transforms into the modulated tilted vortices. The 
modulation vanishes at a second-order phase transition at which the system 
transforms into the tilted-chain state. This phase transition provides possi- 
ble interpretation of the recent Lorentz-microscopy experiment^ at which it 
was observed that the pancake stacks located in the chains smear along the 
chain direction while the pancake stacks outside chains remain well defined. 
This smearing may indicate transition into the tilted-chain state. 

Another phase transition occurs at very small densities of pancake vor- 
tices, a ~ 20 — 30Aj, and only when a exceeds a certain critical value ~ 0.5 
(exact criterion depends on the in-plane magnetic field). In this case a 
small c-axis field penetrates in the form of kinks. The kinked vortex lines 
forming tilted chains are composed of JV pieces separated by kinkP^^. If 
the kink energy is only slightly smaller then the energy per pancake in a 
straight pancake stack then at very small concentration of kinks, typically 
at o ~ [20 — 30] A J, the kinked chains are replaced with strongly deformed 
crossing chains via a first-order phase transition. Due to the opposite signs 
of interactions (kinks repel and deformed pancake stacks attract each other) 
this transition is accompanied by a very large jump in the pancake density. 
With further decrease of the pancake separation the chain smoothly trans- 
forms back to the tilted chain as it was described in the previous paragraph. 

We constructed the chain phase diagrams for different ratios a. As 
follows from the above description, there are two types of phase diagrams in 
the region of intermediate a's: 

• In the range 0.4 ^ a ^ 0.5 a small c-axis field first penetrates in 
the form of pancake-stack chains located on J Vs. Due to attractive 
coupling between the deformed stacks, their density jumps from zero 
to a finite value. With further increase of the c axis field the chain goes 
through intermediate strongly deformed configurations and smoothly 
transforms into a tilted chain; 

• In the range 0.5 ^ a ^ 0.65 small c-axis field first penetrates in the 
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form of kinks creating kinked tilted vortices. With increase of a c axis 
field this structure is replaced via the first-order phase transition with 
the strongly-deformed crossing chains. This transition is accompanied 
by a large jump of pancake density. Further evolution of the chain 
structure is identical to the smaller-a scenario: the structure smoothly 
transforms back into the tilted-chain state. 

The exact transition between the two types of behavior depends also on 
in-plane field. Larger in-plane field favors the first scenario. 

The second scenario provides a natural interpretation to recent scanning- 
Hall-probe observation^^. It was observed that at very small concentration 
of pancakes the chains are magnetically homogeneous and separate pancake 
stacks are not resolved. When the external field exceeds some critical value 
of the order of several Oersted, crystallites of the pancake stacks are sud- 
denly formed along the chain and the flux density in the crystallites is much 
higher then the flux density in the homogeneous chains. Our calculations 
provide consistent interpretation for these observations. The magnetically 
homogeneous chains are interpreted as kinked/tilted chains (such interpre- 
tation has been proposed by Dodgson^'^) and formation of crystallites can be 
attributed to the low-density [kinked lines] -[crossing chains] first-order phase 
transition with density jump. 

The paper is organized as follows. In Section |3 we consider the chain 
energy functional. In SectionEDwe present analytical estimates for the chain 
energy in the two limiting cases, crossing and tilted chain, and discuss lo- 
cation of transitional region between these states. In section |^ we review 
attractive interaction between deformed pancake stacks located on JVs. Sec- 
tion ED contains the results on numerical exploration of the the phase dia- 
gram. We discuss two different phase transitions between the tilted and 
crossing chains and two types of phase diagrams which are realized in the 
region of intermediate parameter a. 



2. CHAIN ENERGY FUNCTIONAL 

In this paper we focus on the structure of an isolated vortex chain with 
period a in x direction and period c = ISl s in z direction where = 2Ni is the 
number of layers per vertical chain period (see Figure^. This corresponds 
to the tilting angle d of the magnetic induction with respect to the c axis 
with V = tan0 = a/c. The pancake separation a is determined by the c-axis 
field Bz- Note that we consider the region of very small (0.1-5 G), where 
it is typically much smaller than external field. The vertical period c is fixed 

by the in-plane field B^, c Jl^^ji^-iB^). For BSCCO (7 ~ 500) this 
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period is approximately equal to 20 layers at 50G. The chains are 

separated by distance Cy, Cy = ^q/cBx ^ c. 

Our calculations are based on the Lawrence-Doniach free-energy func- 
tional in the London approximation, -FlDj (see, e.g., Ref. The chain 
structure is mainly determined by pancake displacements from aligned po- 
sitions Un and regular phase distribution (pr,n{x,y)- We consider the case 
c <C A and in-plane distances much smaller than Ac- In this situation the 
general Lawrence-Doniach energy, F^u, can be significantly simplified us- 
ing several approximations: (i)we neglect screening of regular phase and z- 
axis vector-potential; (ii)we subtract the energy of straight pancake stacks, 
(Sa;/$o)^PS) allowing us to eliminate logarithmically diverging pancake-core 
contributions; (iii)we drop the trivial magnetic energy term B'^/Stt which 
plays no role in selection between different chain phases. We will use the 
chain energy per unit area lattice, E = Cy[Fu:)/V — B'^/^Stt) — {Bz/^o)sps] 
with V being the total system volume. With the above assumptions the 
chain energy functional can be represented as 
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+ 7;T7~yZ^Mriun-Um,n-m) (1) 

where the phase stiffness, J, and the Josephson coupling energy, Ej, are 
given by 

J = — , Ej = 2' with eo = 7^^, (2) 

A = Xab and Ac are the components of the London penetration depth and 7 = 
^c/^ab- The ratio of the two energy scales determines the most important 
length scale of the problem, the Josephson length, Aj = 7s = a/ J/ Ej . 
(j)v,n (r) is the vortex phase variation induced by displacement of pancake 
rows, Un, from the ideally aligned positions 

tan {tt{x + Un)/ a) tan(7rx/a) 

ipv,n{x, y; Un) = aictan — - — — arctan- — — --; (3) 

tanh(7ry/aj tanh{TTy/a) 

the discrete gradient S/n4>n is defined as Vn4>n = (pn+i — 0n- UMr{un — 
Um,n — ITT') is the interaction energy between the pancake rows per unit 
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length, computed with respect to straight stacks 

UMr{u, n) = - ^ [C/A/(ma + u,n) — UMima, n)] 



where C/m(R, n) is the magnetic interaction between pancaked 
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Ei{u) = (exp(— dv is the integral exponent, r = ^JE? + {ns^ , and 
L is a cutoff length; Ntot is the total number of layers. 

The energy contains contribution coming from a long-range suppres- 
sion of the Josephson energy accumulated from distances c <C y <C Cy, that 
is identical in all chain phases, -Ej-lr = Ej7r'^Cy/{6N'^s). It is convenient to 
separate this term and define the local chain energy, E''°'^, 
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(5) 



This part of energy weakly depends on Cy and does not diverge for Cy ^ oo. 

We numerically minimized the energy functional with respect to 
and (pr,n{x-,y) within unit cell, 0<x<a, 0<y< Cy/2, and < n < 
Ni = N/2 using appropriate symmetry and periodic boundary conditions. 
In calculation of magnetic coupling energy one has to take into account 
periodic conditions for pancake displacements, Un+N = Un- In addition, if 
one selects z axis origin at the center of the JV then symmetry also requires 
U-n = —Un- Two simple limiting cases in Fig. ^correspond to (i) Un <^ a 
for the crossing chain and (ii) n„ = — a(l — (n — l/2)/Ni)/2 for the tilted 
chain. 



3. CHAIN ENERGIES IN LIMITING CASES 

Analytic estimates for energy contributions are possible in two limiting 
cases of weakly deformed crossing chain and tilted chain shown at Fig. 
In this section we summarize results for chain energies in different limits. 
Detailed derivations will be published elsewhere. 

Local energy of crossing chain per unit area is given by 



(6) 
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Here the first term represents the local JV energy, i^jy , and the second term 
is the contribution from the crossing energies of JVs and pancake stack^. 

The local tilted-chain energy per unit area in the limit v = tan^ <C 7 
is given by 



TTitl 



a \ 



(7) 



with ~ 0.5772 being the Euler constant. In Eq. ((71) the first term rep- 
resents the loss of the magnetic coupling energy in the tilted chain and the 
second term represents the Josephson energy loss. 

In the region of kinked lines ^ 7 (but v <^ Nj/2tt), the local chain 
energy is given by 



rploc ^ rploc _i_ 



0.44 ^ / /7^\ 



ln-- + -^ln^-J- 2.454 



(8) 



In this equation -Ejy is the local JV-lattice energy, the first term in square 
brackets gives the single-kink energy and the third term gives the kink in- 
teraction energy. 

To find out whether the crossing or tilted chain is realized for given 
values of parameters a, c, and a, we have to compare the energies of these 
states. Comparison of Eqs. © and gives transitional pancake separation 
Oc which decreases with and increases with a Naively, one may think 
that intersection of the energy curves for two states would correspond to a 
first-order phase transition between these states. However, as we will see 
from numerical simulations, in the region 1^ = tan 9 < 7 another scenario is 
realized. Typically, strongly deformed intermediate chain configurations de- 
velop in the transitional region providing a smooth transformation between 
the two limiting configurations. Therefore, a simple energy comparison gives 
only an approximate location of the transitional region separating the two 
configurations. 

In the region 1 < Z//7 < N/2ti and a > 2tt\ comparison of the energies 
(jSj) and (jHl) gives the following equation 

This equation has a solution only in the kink penetration regime. In (1/a) — 
0.81 < 0, where the kink energy is only slightly smaller that the energy per 
pancake of a straight pancake-vortex stack. In contrast to the case 1/ <C 7, 
this equation does correspond to a very strong first-order phase transition. 
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4. ATTRACTION BETWEEN DEFORMED PANCAKE 
STACKS. MAXIMUM EQUILIBRIUM SEPARATION 

A peculiar property of the crossing chain is an attractive interaction 
between the deformed pancake stacks at large distance^^. As a consequence, 
when the magnetic field is tilted from the direction of layers the density of 
pancake stacks located on JVs jumps from zero to a finite value. This means 
the existence of a maximum equilibrium separation am between pancake 
stacks, i.e., chains with a > am are not realized in equilibrium. Note that 
the tilted vortices also attract each other within some range of angles and 
distances ^ meaning that tilted chains also have this property in some range 
of parameters. 

Simple analytical formula for the attraction energy between the de- 
formed pancake stacks can be derived for very anisotropic superconduc- 
tors Aj ^ A in the range A ^ i? ^ Aj In this limit short-range 
pancake displacements Un from the aligned positions in two neighboring 
stacks produce dipole-like contribution to interaction energy per unit length 
between these stacks 5Ui{R) = — 2eo ("u^) / ■ This term has to be com- 
bined with the usual repulsive interaction between straight stacks Uio{R) = 
2eoKo(-R/A) « 2eo\/7rA/2i?exp (— i?/A). Minimum of the total interaction 
energy, Uio{R) + 6Ui{R) gives an estimate for the maximum equilibrium sep- 
aration am^, dm = Aln(CA^/ {u"^)) and this result is valid until am < Aj. 
Because in BSCCO Aj is only 2 — 3 times larger than A, this simple for- 
mula is not practical for this compound. We will see that am in BSCCO is 
usually larger than Aj. In general, the maximum equilibrium separation am 
is determined by the minimum of the pancake part of energy per one stack 
U{a) = a{E{a) — E{oo)). When the main contribution to the total inter- 
action energy is coming from the nearest -neighbor interciction, coincides 
with the position of the minimum in the pair interaction potential. 



5. NUMERICAL EXPLORATION OF CHAIN STRUCTURES 



We explored the chain structures by numerically minimizing with 
respect to the pancake row displacement n„ and regular phase distribution 
(pr,n for different values of the parameters a, N = 2Ni, and a. In the following 
sections we review the results of these calculations. 
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5.1. Phase Transition Prom Crossing To Tilted Chains With 
Decreasing Pancake Separation 

We studied the evolution of chain structures with decreasing pancake 
separation a at fixed a and A^. For small values of a we found that the 
chain structure smoothly evolves with decreasing a from crossing to tilted 
configuration. An example of such evolution is presented in Fig.|21for = 14 
and a = 0.4. The upper plot shows the dependence of the maximum pancake 
displacement from the straight-stack position Umax/cL (defined in the inset) 
on the pancake separation a. We will use this parameter to characterize 
the chain structure throughout the paper. It changes from zero for straight 
stacks to (1— l/iV)/2 for tilted chains. The lower plot shows the a dependence 
of the reduced pancake energy per unit cell, U = ca{E{a) — E{oo))/ J, which 
is used to determine the maximum equilibrium separation Um- At large a a 
weakly deformed crossing configuration is realized (see structure at a = 3Aj). 
With decrease of a the chain evolves into strongly corrugated configurations 
such as configuration for a = 2Aj in Fig. [21 With further decrease of a 
this structure smoothly transforms into modulated tilted lines (see structure 
for a = 1.5Aj). Finally, the last structure transforms via a second-order 
phase transition into the straight tilted lines. For parameters used in Fig. |21 
this occurs at a = 1.3Aj. The plateau in the dependence Umax{a)/a below 
this value of a corresponds the maximum possible relative displacement (1 — 
l/N)/2 w 0.4643 in the tihed chain. 

Comparison of numerical and analytical calculations shows that the 
analytical estimates © and ((TJ accurately reproduce numerical results for 
the weakly deformed crossing chain and for the tilted chain. However, the 
numerical study predicts intermediate configurations with energies smaller 
than the energies of the both limiting configurations. As a result, a naively 
expected first-order phase transition is replaced by a continuous transition 
occurring at significantly smaller a. The transitional region just marks the 
location of the intermediate strongly deformed chain configurations. The 
observed continuous phase transition is relat ed to the instability of isolated 
tilted vortices in anisotropic superconductor j ^^ l ^^ l 

In Fig. 1^1 we present the chain phase diagrams in the a-A^-plane for 
a = 0.4. The thick line shows the phase transition into the tilted-chain 
state. One can see that at larger the transition takes place at smaller 
a. With increase of a this line moves higher. Dotted lines show locations 
of the maxiiimm eQuilibrium scpctration (ifji discussed in Sec. El We see 
that am{N) line crosses the transition line meaning that at small N falls 
into the tilted-chain region and at large it falls into the crossing-chain 
region. The obtained phase diagram implies that at small tilting angle of 
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Fig. 2. Upper plot shows the dependence of the maximum displace- 
ment (defined in the inset) on pancake separation a for = 14 and a = 0.4. 
Lower plot shows dependence of the pancake energy per unit cell U (in units 
of phase stiffness J) on pancake separation. Its minimum determines the 
maximum equilibrium separation am- Representative chain structures are 
shown for several values of a(boxes show positions of the pancake vortices 
and horizontal lines mark JVs). One can see that the system evolves from 
weakly deformed chain (a = 3Aj) via strongly deformed chain (a = 2Aj) 
to modulated tilted chain (a = 1.5Aj). The last structure transforms via a 
2nd-order phase transition at q = 1.3Aj into tilted straight vortices. 



a = 0.4 



—3 




°8 10 12 14 16 18 

N 



Fig. 3. Chain phase diagrams in the plane a/Xj — N for the ratio a = 0.4. 
Thick line with open circles indicates phase transition into the tilted chain. 
Thin line shows the nicixirnuin Gqiiilibriiiiii sGpara-tion drw 
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the field (corresponding to small a) the tilted chains have lower energy than 
the crossing chains. This is similar to the situation at higher fields, in the 
dense lattice, where the crossing-lattices state also is expected to transform 
into the simple tilted lattice at small tilting angle of the fielcP. 



5.2. Reentrant Transition To Kinked/Tilted Lines At Small 
Concentration Of Pancakes At a > 0.5 

At higher values of the ratio a a new qualitative feature appears in 
the phase diagram. When a exceeds the critical value, the c axis field ini- 
tially penetrates in the form of kinks forming kinked vortex lines (lock-in 
transitiorPEI^ . The critical value of a is determined by combination of nu- 
merical constants in the pancake-stack and kink energies. At present, our 
best estimate for this constant is Uc ~ 0.44. The critical value of a in- 
creases with decrease of N, due to the increasing relative contribution of the 
crossing energies in the energy of crossing chain. A very peculiar behavior 
is expected when a only slightly exceeds ac- The competing chain states 
have very different interactions: deformed stacks attract and kinks repel 
each other. Moreover, at the same value of the c-axis magnetic induction, 
Bz , the kinks are much closer than the stacks and the absolute value of the 
kink interaction energy is much larger than the interaction energy between 
deformed stacks. As a consequence, with increase of the total energy of 
the kinked lines rapidly exceeds the total energy of the crossing chain and 
the system experiences a first-order phase transition into the crossing-chain 
state. Due to the attractive interaction between the pancake stacks, the 
pancake/kink separation at which the energy curves cross does not give the 
equilibrium separation for the crossing chain and the stack separation jumps 
at the transition to a value slightly smaller than the maximum equilibrium 
separation am- This means that the phase transition is accompanied by 
jumps of the pancake density and magnetic induction, B^. 

This behavior was confirmed by numerical calculations. Figure |1] shows 
a plot of the dependence of the energy per chain unit cell U on the pancake 
density n = 1/a for A = 16 and a = 0.6. This dependence has two branches, 
corresponding to two different starting states at small n, crossing chain and 
kinked lines, and the kinked vortex lines have smaller energy at smaller n. 
The branches cross at n = 0.068 marking the first-order transition. The 
variations of U at small n occur due to the interaction energy and one can 
see that the kink interaction energy is much larger than the crossing-chain 
interaction energy. With further increase of n, the crossing chain smoothly 
transforms into the tilted chain following the scenario described in the pre- 
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Fig. 4. The dependence of the pancake energy per chain unit cell on pancake 
density n = 1/a for = 16 and a = 0.6. The two branches at low n 
correspond to two different starting states, kinked lines (1) and crossing 
chain (2). The kinked lines have lower energy at very small n, at n < 0.068. 
Crossing chain smoothly transforms back into tilted chain with increase of n. 
The transformation is completed at the second-order phase transition point 
near nXj = 0.4. Chain configurations are shown at four marked points. 

vious Section. The second-order phase transition for these parameters takes 
place at n ~ 0.4 corresponding to a ~ 2.5Aj. The latter value is somewhat 
smaller than the maximum equilibrium separation ~ 3.44Aj. 

The pancake density in the chain can not be directly fixed in experiment. 
Instead, the magnetic field strength, H^, fixes the chemical potential fin, 
Hh = ^o-f^z/(4vr), and the equilibrium density is determined by the global 
minimum of the thermodynamic potential G{n) = E{n) — ^nn. To find 
the density evolution with increasing the chemical potential, we plot in Fig. 
Efleft) the density dependencies of the reduced thermodynamic potential, 
5£—^n for different /i and representative parameters = 14 and a = 0.6. As 
the energy of isolated stacks is subtracted in the dimensionless chemical 
potential is shifted with respect to its bare value and it is related to the 
magnetic field strength as /x = N^o{Hz — ifci)/(4eo), where Hd is the lower 
critical field for H||c. We find that for selected parameters the transition 
takes place aX ^ = = —6.668. At /i < /i^ the global minimum falls into 
the region of kinked lines and at /x > /xj it jumps into the region of crossing 
chain. Note that the density value, at which the energy curves cross (kinks 
in the lines in Fig. , is always larger than the lower density from which the 
jump to the high-density state takes place. At the transition, the density 
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Fig. 5. Left plot shows the density dependence of the pancake part of the 
thermodynamic potential per unit area, 6£ — jin (in units of eo/(vrA jN)) at 
different values of chemical potential fi corresponding to different external 
fields for = 14 and a = 0.6. Kinks in the curves separate regions of 
tilted/kinked lines (low n) and crossing chains (high re). The equilibrium 
density is given by the global minimum of this energy. One can see that at 
fi ~ —6.668 the system experiences a first-order phase transition with very 
large density jump. Right plot shows the corresponding dependencies of the 
average magnetic induction, Bz, (left axis) and the maximum induction in 
the middle of the chain, Bch,z, (right axis) vs the magnetic field strength, 
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a = 0.6 




N 

Fig. 6. The left panel shows the phase diagram in the N-a plane for a = 0.6. 
The right panel illustrates definitions of the phase lines using the plot of the 
thermodynamic potential 6£ — vs a/Aj at the transition point. The lines 
1 and 1' correspond to the two limiting pancake separations at the transition 
point between which the jump occurs. The line 2 indicates crossing of the 
energy curves for the kinked and crossing chain. The line 2' shows position 
of the continuous transition into the tilted chain. Dotted line slightly above 
I'-line shows position of the maximuiTL 6Qiiilibriiim separation cirn' 

We also 

show by the dashed line the crossover line a/Xj = N above which well- 
defined kinks appear. 
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jumps almost ten times, from 0.037/Aj to 0.315/Aj. The right plot in Fig. 
[5] shows the corresponding dependencies of the average magnetic induction, 
Bz = ^on/cy, (left axis) and the maximum induction in the middle of the 
chain, Bch,z = ^on'/^i (right axis) vs the magnetic field strength, Hz- One 
can see that the jump in Bz occurs from ~ O.IG to ~ IG. 

The numerically obtained phase diagram in the N-a plane for a = 0.6 is 
shown in the left panel of Fig. El The plot of the thermodynamic potential 
at the transition point in the left panel illustrates definitions of different 
lines in the phase diagram. The lines 1 and 1' show the limiting pancake 
separations at the first-order transition between which the jump takes place. 
When the chemical potential is fixed by external conditions the area between 
these lines is bypassed in equilibrium. The line 2 inside this region marks the 
crossing of the energy curves for the two states. At large the jump takes 
place from the kinked- lines state into the strongly corrugated configuration. 
This configuration transforms into the tilted chain with further decrease of 
a via continuous transition shown by the line 2'. Below = 14 only tilted 
chains realize, but the density jump still exists. The upper separation grows 
approximately proportional to A^ while the lower separation slowly decreases 
with A^ and lies slightly below the maximum eQuilibrium separation ci^ 
shown by a dotted line. This means that the relative density jump increases 
with A^. This type of phase diagram exists within a finite range of a, roughly 
0-5 ^ Oi < 0.7. At larger a's only tilted chains exist. 



6. CONCLUSIONS 



In conclusion, we investigated the phase diagram of an isolated vortex 
chain in layered superconductors. In the region where Josephson and mag- 
netic coupling are approximately equal, we found a very rich behavior. The 
crossing chains typically transform into the tilted chains with decreasing 
pancake separation via formation of intermediate strongly deformed config- 
urations and continuous phase transition. When the relative strength of the 
Josephson coupling exceeds some typical value, the phase diagram becomes 
reentrant. At very small c-axis field, tilted chains are realized in which the 
vortex lines have kinked structure. With increasing the c-axis field these 
low-density tilted chains transform via a first-order phase transition into the 
strongly-deformed crossing chains. This transition is accompanied by a large 
jump of pancake- vortex density. With further field increase these crossing 
chains transform back into the tilted chains via a second-order transition. 
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